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Abstract. We report the detection of a cyclotron absorp- 
tion line and its second harmonic in the average spectrum 
of the high mass X-ray binary 4U1907-f09 observed by 
the BeppoSAX satellite on 1997 September 27 and 28. 
The broad band spectral capability of BeppoSAX allowed 
a good determination of the continuum against which the 
two absorption features are evident at ~ 19 and ~ 39 keV. 
Correcting for the gravitational redshift of a ~ 1.4 
neutron star, the inferred surface magnetic field strength 
is Bsurf = 2.1 X 10^^ G that is typical for this kind of sys- 
tems. We also detected an iron emission line at ~ 6.4 keV, 
with an equivalent width of ^ 60 eV. 

Key words: binaries: close — star: neutron — pulsars: 
individual: 4U 1907+09 — X-rays:stars 



1. Introduction 

4UI907-I-09 is an X-ray pulsar powered by wind accretion 
from a close blue supergiant companion star. It belongs to 
the class of the high-mass X-ray binaries (see, e.g., White, 
Nagase & Parmar 1995). 

4UI907-I-09 was discovered in the Uhuru surveys (Giac- 
coni et al. I97I, Schwartz et al. 1972). An observation of 
this source performed by the RXTE satellite has been re- 
cently analyzed by in 't Zand et al. (1998). The reported 
pulse period is 440.341 ±0.006 s. No pulsations were found 
above 40 keV up to date. The reported orbital period is 
8.3753 ± 0.0001 d and the eccentricity is 0.28 ± 0.04. The 
pulsar showed a fairly constant spin-down between 1983 
and 1996 with a mean derivative of P = -1-0.225 s yr~^. 
Periodic flares, two per orbit, have been reported (Mar- 
shall & Ricketts 1980) with an increase of the intensity 
by about a factor 5. The energy spectrum of 4U1907+09 
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was observed with several X-ray missions (Ariel V, Tenma, 
EXOSAT, Ginga, XMPC and RXTE) covering, not simul- 
taneously, the range 2-100 keV. It was fitted by an ab- 
sorbed power law (Nh in the range 1.5 — 5.7 x 10^^ cm~^ 
and photon index in the range 0.83 - 1.52 ~ see Schartz 
et al. 1980, Marshall and Ricketts 1980, Makishima et al. 
1984, Cook and Page 1987, Chitnis et al. 1993). Data from 
a Ginga observation of 4U1907+09 were fitted using the 
NPEX model (Negative and Positive power-laws Expo- 
nential) as continuum plus a Lorentzian shaped cyclotron 
absorption line at ~ 19 keV and an iron line at 6.6 keV 
with a flux of 1.7±0.5 ph s"^ In the RXTE observation 
the spectrum of 'non dipping' data, subtracted by the 'dip- 
ping' spectrum, was fitted in the 2-30 keV range by an 
absorbed power law with high energy cut-off at 13.6 keV 
with no evidence of any absorption features (in 't Zand et 
al. 1997). 

In this paper we report the results of a spectral analy- 
sis of the broad band spectrum (1-80 keV) of 4U1907-H09 
observed with the BeppoSAX Narrow Field Instruments 
(NFIs) . We confirm the presence of the absorption feature 
at ^ 19 keV and we report the detection of a strong second 
harmonic. 

2. Observation 

4U1907-I-09 was observed by the NFIs on board Bep- 
poSAX (Boella et al. 1997a) on 1997 September 27 and 
28. The NFIs consist of four coaligned instruments: the 
Low Energy Concentrator Spectrometer (LEGS) operat- 
ing in the energy range 0.1-10 keV (Parmar et al. 1997), 
the Medium Energy Concentrator Spectrometer (MEGS) 
having three units operating in the range 1-10 keV 
(Boella et al. 1997b), the High Pressure Gas Scintillation 
Proportional Counter (HPGSPC) operating in the range 
4-120 keV (Manzo et al. 1997) and the Phoswich Detec- 
tor System (PDS) with four scintillation units operating 
in the range 15-300 keV (Frontera et al. 1997). The en- 
ergy resolution (AE/E) is 8% at 6 keV, 6% at 20 keV and 
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12% at 40.0 keV. We do not consider LEGS data because 

the MEGS statistics is more than enough to model the 
continuum under the cyclotron lines. The total exposure 
time was 68576 s for the MEGS, 33079 s for the HPGSPG 
and 34731 s for the PDS. HPGSPG and PDS were exposed 
about half as long as the MEGS due to the rocking colli- 
mator observation mode. The MEGS events were reduced 
using the SAXDAS v. 1.1.0 package (BcppoSAX Gook- 
book, http://www.sdc.asi.it/software/cookbook) while 
HPGSPG and PDS data were reduced using the XAS 
v. 2. 0.1. package (Ghiappetti & Dal Fiume 1997). Response 
matrices were used which were made publicly available on 
August 31, 1997. 

The 4U1907+09 background subtracted light curves 
of the three NFIs are reported in Figure 1. The gaps in 
the light curves are due to non-observing intervals dur- 
ing passes over the South Atlantic Anomaly and during 
Earth occultations. An increase in intensity of about a 
factor two (up to 13 c/s, ~ 45 mGrab in the MEGS) is 
observed at ^ 2500 s after the beginning of the observa- 
tion. Using the Porb reported by in 't Zand et al. (1998), 
the flux enhancement observed by BeppoSAX is found to 
occur 613.93 ±0.05 orbital periods after the primary flare 
observed by Tenma (Makishima et al. 1984). Therefore we 
tentatively identifled this episode with the periodic pri- 
mary flare, confirming the occurrence of the phase-locked 
flares. 



MECS 
1.6— lO keV 




Time from start (sec) 



Fig. 1. Light-curves of the NFIs binned at the spin period of 
4U1907-I-09 



3. Spectral Analysis 

Spectral analysis was performed on the spectrum aver- 
aged over the complete observation. Due to the presence 



of the galactic ridge and the supernova remnant W49B 
0.8° from 4U1907+09) additional care has been ap- 
plied in the determination of the background. The local 
background has been measured in a region of the image 
20 arcminutes away from the source for the MEGS, and 
compared with the background of the archive blank fields. 
The MEGS local background spectrum can be well repre- 
sented by the blank field increased by a factor 1.6. There- 
fore blank sky background has been used for the MEGS, 
applying the above correction. For the coUimated instru- 
ments we use the local background accumulated in the 
off-source position. Anyway, from previous measurements 
of W49B (Fujimoto et al. 1995, Smith et al. 1985) and the 
galactic ridge (Yamasaki et al. 1997) and considering the 
strong reduction in the effective area for off axis sources 
we estimate that the contribution of these sources is neg- 
ligible above 10 keV. 

The energy range used in the spectral analysis for each 
NFI was: 1.6-10 keV for MEGS, 10-30 keV for HPGSPG 
and 15-80 keV for PDS. The continuum in the range 1-80 
keV has been modeled adopting an absorbed power law 
with an exponential cut-off for energies above Ecut- To 
take into account the systematic differences in the normal- 
ization between the different NFIs (Gusumano et al. 1998) 
the relative normalizations have been kept as free parame- 
ters. The Xred value with the above continuum model was 
2.14 (221 dof). Figure 2 shows the spectra with the best 
fit continuum and the residuals. 




2 5 10 30 SO 

energy (keV) 



Fig. 2. 1.6-80 keV spectrum of 4U1907+09 with the best fit 
continuum model (top) and residuals (separately for each NFI 
in the bottom three panels). 

The excess at ~ 6.4 keV in the residuals strongly sug- 
gests the presence of an iron emission line. Therefore we 
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servation was 4 x lO^^D^^^ in the 2-10 keV energy range 

in the 2-80 keV range where Df^^, the dis- 
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Fig. 3. 1.6-80 keV spectrum of 4U1907+09 with the best fit 

continuum model plus an emission gaussian line at --^ 6.4 keV 
and two absorption gaussian lines at 18.8 and 39.4 keV. The 
residuals do not show significant features. 



included an additional gaussian to model the iron line. 
Other features are evident from the residuals above 10 
keV. A feature at ^ 40 keV is evident. We tried several 
models to fit this feature and we obtained the best fit 
including an absorption according to a gaussian shaped 
feature (Soong et al. 1990) centered at ^ 39 keV. The 
Xred value after the addition of the absorption gaussian 
feature was 1.27 (215 dof). A feature at ~ 19 keV is also 
present in the residuals of Fig. 2. We tried to improve 
the fit adding an absorption gaussian factor to the above 
model (model 1). The Xred value after addition of this 
last gaussian shaped feature was 1.06 (212 dof). Figure 
3 shows the result of the fit. The residuals do not show 
significant structures. The F-test for three additional pa- 
rameters gives a significance of the absorption feature at 
~19 keV greater than 99.999%. The significance of the last 
feature and its energy that is compatible (Table 1) with 
half of the energy of the other absorption feature allow 
us to identify these two absorption features as the funda- 
mental and the second harmonic of the cyclotron absorp- 
tion line at ~ 19 keV reported by Makishima and Mihara 
(1992) from Ginga observations. The Lorentzian absorp- 
tion model of Mihara (1995) (model 2) was also tried to 
fit the two absorption features. Also in this ease the fit 
gives an acceptable x^- Best fit parameters obtained with 
the models described above are reported in table 1. 
Other descriptions of the continuum (such as NPEX, Mi- 
hara 1995), together with the models for cyclotron lines 
described above, gave worse results. 
The luminosity of the source during the BeppoSAX ob- 



tance of 4U1907-I-09 in kpc, has been roughly estimated 
between 2.4 and 5.9 kpc (Van Kerkwijk et al. 1989). 



Table 1. best fit parameters 
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Ecut (keV) 


12.0 ± 0.3 


12.3 ±0.4 


E/oid (keV) 


12.0 ± 0.3 


13.8±0.8 


EcflFi (koV) 


18.8 ± 0.4 


19.3±0.2 


TCijFi (keV) 


2.2 ± 0.4 




EWcRFi (keV) 


2.3 ± 0.4 




DcRFl 




0.29±0.05 


WcRFi (keV) 




1.8±0.8 


EcHF2 (keV) 


39.4 ± 0.6 


2 X EcflFi 


CTCRF2 (keV) 


3.6 ± 0.7 




EWcflF2 (keV) 


16.7 ± 2.1 




DcRF'2 




2.7±1.1 


WcRF2 (keV) 




2.8±1.3 


Xrid (d.O.f.) 


1.06 (212) 


1.03(213) 



Norm is in unit of photons/keV/cm^/s a 1 keV 
N_H is in unit of (10^^ atoms cm~^). 
I is the total photons cm^^ s^^ in the line. 
Quoted errors refer to single-parameter 68% confidence level 
The upper limit to the Oiron is at 90% confidence level. 



4. Discussion 

The detection of cyclotron lines in the hard X-ray spec- 
trum of binary pulsars allows a direct measurement of the 
magnetic field of strong magnetized neutron stars. Ginga, 
with its relatively broad energy band (1-40 keV), allowed 
the detection of cyclotron absorption lines in several X-ray 
pulsars. Sometimes also the presence of a second harmonic 
was claimed. However the detection of these harmonics 
was often ambiguous, as in the case of 4U1538-52 (Clark 
et al. 1989). In fact only the red wing of the second har- 
monic (at ~ 40 keV) could be seen with Ginga and a high 
energy cut-off could also fit the spectrum. The very broad 
energy band of BeppoSAX NFIs allows to measure the X- 
ray spectrum of 4U1907±09 from 1 to 80 keV, well above 
the top end energy of Ginga Large Area Counter. So we 
can constrain the shape of the underlying continuum with 
great accuracy, allowing an easier measure of spectral fea- 
tures, when present. In the case of this source, we detect 
two absorption features at 18.8 keV and 39.4 keV respec- 
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tively, which we interpret as two cyclotron lines, the fun- 
damental and the second harmonic. In this case one might 
expect that the line at 19 keV would be deeper than the 
line at 39 keV because the cyclotron opacity decreases for 
higher harmonics. On the contrary for this source the sec- 
ond line appears deeper than the first line. A similar line 
structure was found in the spectrum of GB 880205 (e. g. 
Alexander and Mcszaros 1989). According to Alexander 
and Meszaros (1989) this line structure can be explained 
through the process of Compton scattering with multiple 
photon emission. In fact, in the multiple emission case, 
each scattering removes a photon with E ~ 2Ec and gives 
two photons with E Ec- Thus photons at ~ 2Ec are de- 
pleted, while those at ^ Ec replenished, and the line at the 
second harmonic becomes deeper than the fundamental. 

Considering the relation between the cyclotron energy 
and the magnetic field Ec/{1 keV) = 11.6 B/{10^'^ Gauss), 
the observed value of the cyclotron energy of 18.8 keV im- 
plies a magnetic field of Bobs — 1-6 x 10^^ Gauss. If the 
cyclotron absorption takes place near the neutron star sur- 
face, where the magnetic field is strong, the observed reso- 
nance energy will be affected by the gravitational redshift: 
Ef^ =Ec{l + z)-'^, with 

where M^s is the mass of the neutron star and R is the 
distance of the region where the line is formed from the 
center of the neutron star. Using Alys = 1.4M0 and R 
111" cm, we get (1 + z)-'^ = 0.76 and E^ ~ 24.7 keV. In 
this case the magnetic field should be Bsurf — 2.1 x 10^^ 
Gauss. 

An estimate of the temperature of the region in which 
the cyclotron absorption is formed can be obtained from 
the width of the observed lines. In fact if the electrons have 
a temperature kT, the energy of the scattered photons 
distributes in a gaussian-like peak with Lorentzian wings 
(Voigt function). In this case the doppler broadening of a 
line /S.E/E is constant and is given by (see e.g. Rybicki 
and Lightman 1979): 

AE = E^J^ (2) 

a factor cos^ - where 9 is the angle between the mag- 
netic field lines and the line of sight - arises if the plasma 
motion is confined by the field. Averaging this over angles 
between and it/ 2 introduces corrections ~ 1. We observe 
that the obtained width of the second harmonic is double 
that of the fundamental, in agreement with the relation 
above. The obtained values of the parameters give a tem- 
perature of 4 - 7 keV. 
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